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Abstract. Leptogenesis is studied within the seesaw neutrino mass model in a regime 
where all sterile neutrinos have prompt rather than delayed decays. It is shown that 
during neutrino thermal production lepton asymmetries are generated in both active 
lepton and sterile neutrino sectors. The large B — L asymmetry is slowly pumped into 
the chemically decoupled right-handed quarks and leptons and baryon number sector 
which later protect B — L from fast L violating processes. The dependence of the final 
baryon asymmetry on couplings and masses is totally different from the decay scenario. 
B does not vanish in the limit of degenerate light neutrinos and the observed asymmetry 
is naturally obtained for a sum of square masses, m 2 , between the atmospheric neutrino 
mass gap and 0.2 eV 2 . 

1 Introduction 

Leptogenesis is an attractive way of generating the baryon number of the Uni- 
verse JI]. The main idea put forward by Fukugita and Yanagida [2] is that a lep- 
ton number asymmetry can be produced in the decays of heavy sterile neutrinos 
into leptons and Higgs bosons and such an asymmetry is partially transferred 
to the baryon sector through electroweak sphaleron processes [3] that violate 
B and L but not B — L. The mechanism requires lepton number and CP vi- 
olation, ensured by neutrino Majorana masses and complex Yukawa couplings. 
Both masses and couplings form the well known seesaw model 0] of light neu- 
trino masses and thus establish a close relation between baryogenesis and low 
energy phenomenology. 

The third Sakharov condition [S] , the departure from thermal equilibrium, is 
satisfied in the delayed decay scenario. It means here that the Majorana neutrinos 
live longer than the age of the Universe when they approach non-relativistic 
temperatures and their production is Boltzmann suppressed. It turns out that 
such condition puts an upper bound on the natural order of magnitude of the 
light neutrino masses |6llj well below the mass scale indicated by the atmospheric 
neutrino anomaly 7 . This natural bound can be circumvented to some extent 
by fine-tuning the parameters to maximize the produced CP asymmetry but 
even so a quasi-degenerate light neutrino spectrum is excluded |8I9| . 

As the seesaw provides a simple and elegant neutrino mass mechanism, it 
would be nice to have a naturally successful leptogenesis mechanism in the case 
where all heavy neutrinos have prompt rather than delayed decays. Given that 
the heavy neutrino Majorana masses and Yukawa couplings provide the neces- 
sary B—L and CP violations the main concern is to satisfy the out-of-equilibrium 
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condition. This is naturally achieved by assuming that the heavy sterile neutri- 
nos do not couple to the inflaton field and are only thermally produced during 
the radiation era from standard lepton and Higgs particles. Before the heavy 
neutrinos reach thermal equilibrium abundances a non-zero B — L asymmetry 
can in principle be generated. That is the subject of this paper. 

In the next section we describe the mechanism in more detail and in the 
following ones we present the main results concerning first, the generation of the 
primary B — L asymmetry and second, the protection mechanism that prevents 
B — L from being completely washed out when the lepton number violating reac- 
tions reach thermal equilibrium. In the last section we summarize the mechanism 
and make a brief discussion. 



2 Thermal production - prompt decay scenario 

The seesaw mechanism adds to the standard model singlet (left-handed) neu- 
trinos, N a , with heavy Majorana masses and Yukawa couplings to the standard 
lepton and Higgs doublets, k, <p, of the form 

h ia kN a cf> + ^M a N a N a + H.C. . (1) 

Spontaneous breaking of SU(2) x U(l) yields the light neutrino mass matrix 

mij = -(hM- 1 h T ) ij v 2 . (2) 

The proper decay rate of N a into leptons and Higgs is r a = (ftJ h) aa M a /8TT. 
Delayed decay occurs when the ratio to the Hubble expansion rate, K a = r a /H, 
is small at the temperature T = M a . In the radiation era, H = 1.66 gl /2 T 2 /M P , 
where denotes the number of relativistic degrees of freedom, 107.5 in the 
standard model. It is enough to compare the sum 

K = l^K a = 10 _ 3eV (3) 

with the light neutrino scale tr[m] to conclude that the delayed decay condition, 
K a < 1, is in general in conflict with the atmospheric neutrino mass gap [7], 
Am 2 w 3 • 10~ 3 eV 2 , which implies K > 50. Strictly speaking, the delayed 
decay scenario only requires that the lightest of the heavy neutrinos satisfies 
K a < 1, and larger K a values are tolerated for a certain choice of parameters 
that maximize the CP asymmetry |8I9| . But that is not the most natural picture, 
in particular if the light neutrinos are quasi-degenerate |8I9| . 

As it turns out, there is a simple modification of the original leptogenesis 
mechanism capable of realizing the out-of-equilibrium condition and producing 
a large lepton asymmetry. According to the inflation paradigm |10) . the radia- 
tion era starts when a scalar field, so-called inflaton, decays and transforms its 
energy into the form of particles so constituting the initial thermal bath. If ster- 
ile neutrinos do not couple to the inflaton, they can be only thermally produced 
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from standard lepton and Higgs particles through Yukawa interactions. These 
processes are naturally out-of-equilibrium at very high temperatures and it takes 
some time before sterile neutrinos reach thermal equilibrium densities. During 
that period there is a clear departure from equilibrium and a B — L asymmetry 
can be generated. 

In this scenario a primary lepton asymmetry is generated in the lepton dou- 
blet and sterile neutrino sectors. A fraction of it is transferred to the bottom 
6r, tau tr, charm cr and other right-handed fermions through slow Yukawa in- 
teractions. Another fraction is converted into baryon number by weak sphaleron 
processes. It is crucial for the survival of such baryon number and bp, tr, cr, ... 
asymmetries that the heavy Majorana neutrinos decay and vanish from the Uni- 
verse before the weak sphalerons or 6r, tr (cr) Yukawa couplings enter in equi- 
librium, which happens at a temperature around Jl] 10 12 GeV (10 11 GeV). By 
that time the AL — 2 reactions mediated by off-shell Majorana neutrinos should 
also be out-of-equilibrium. In other words, the moment when the AL = 2 reac- 
tions freeze out determines which decoupled sectors or quantum numbers protect 
and store a fraction of the primary B — L asymmetry. Their rate relates to the 
sum of light neutrino square masses rh 2 = tr[m?7v] as 7a ~ fh 2 T 3 / (7r 3 u 4 ) |12I13| . 
They are out of equilibrium, i.e., r 2 /H < 1, at temperatures below 

r ou ,,2^ifY! 10 « Ge v. ,4) 

ni 

We assume of course that the masses M a are larger than T out so that the heavy 
neutrinos vanish above that temperature. T out « 10 12 GeV (10 11 GeV) for neu- 
trino masses m 2 = 0.04 (0.4 eV 2 ), well above the atmospheric neutrino mass 
gap. As will be shown below, if rh 2 < 0.2 eV 2 then the baryon number sector 
and the right-handed quark bp, have the dominant role in protecting B — L and 
producing the present baryon number of the Universe. If m 2 exceeds 0.2 eV 2 
that role would be played by the charm quark cr but not efficiently enough to 
produce the present baryon asymmetry. 



3 Generation of B — L 

We assume that the Universe is initially empty of sterile neutrinos, not produced 
in the inflaton decay but only thermally from standard leptons and Higgs. For 
definitcness we assume a reheating temperature much higher than the neutrino 
masses. The dominant production reaction is the Higgs boson decay into leptons 
and neutrinos, <f> — > liN a , allowed by a Higgs thermal mass, = x^T ~ 
0.6 T |14| . much higher than the lepton thermal masses. The Yukawa top quark 
scattering processes like q{tR — > hN a and other channels contribute to not more 
than 40 % of the total N a production rate. The Higgs boson proper decay rate 
into N a is r£ = (Jv h) aa m^/16TT per isospin state. In the radiation era the 
initially zero neutrino densities n a converge to the fermion equilibrium density 
Ueq = 0.90T 3 /7r 2 with relaxation temperatures T a w 0.35 K a M a : 

n a = n cq (l-e- T ^ . (5) 
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As a consequence of the prompt decay assumption, K a > 70, sterile neutrinos 
reach thermal equilibrium while they are still ultra relativistic. 

Leptogenesis is dominated by the following CP asymmetric reactions: Higgs 
decays into leptons and sterile neutrinos, inverse decays and scatterings of leptons 
off neutrinos. The asymmetries result from the interference between the diagrams 
on the left side of Figure ^ with the absorptive part of the respective diagram 
on the right-hand side. Unitarity and CPT invariance impose constraints on 
the reaction rates such that when all particles are in thermal equilibrium the 
asymmetry source terms cancel each other [15116117] . Generically denoting the 
difference between the total rates of a reaction X — > Y and its CP-conjugate 
X — » Y as Aj(X — > Y) (in a fixed comoving volume), the relevant constraints 
are expressed in terms of thermal equilibrium rates as: 

Ay(0 N a k) eq + Aj(N a k -> 0) cq = , (6a) 
Aj{N a h -► N b l 3 ) cq + Aj(N b l 3 -► NJ t ) cq = , (6b) 

Aj(N a k -» 4>) cq + A~/(N a k -> N b lj) eq = . (6c) 

bj 

The evolution of the set of observables Q a is governed by the Boltzmann 
equations, whose general structure is, integrating over a constant comoving vol- 
ume, Q a = —r a f}Qp + (Qa)s- The first are transport terms and {Q a )s stand 
for source terms. The particle asymmetry sources vanish in thermal equilibrium 
but, as long as the N a neutrinos stay rarefied, scatterings and inverse decays do 
not match Higgs decays and particle asymmetries develop in the various lepton 
flavors, sterile neutrinos, and Higgs boson as well, as enforced by hypercharge 
conservation. The calculation of the source terms is simplified by employing 
Boltzmann rather than quantum statistics as is common practice. Relating the 
out-of-equilibrium and thermal equilibrium reaction rates as 7 = 7 cq n. Q /n oq for 




Fig. 1. Diagrams contributing to the CP-asymmetries in decays, inverse decays and 
scatterings 
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N a X — > Y, and 7 = 7 oq for X — > 7V a y, one derives from the above CPT and 
unitarity constraints the k lepton number sources, 



(U) S = E q - Weq , (7) 

Ay^- - JVftZOeq = -^f^ , (8) 

where c; « 0.1, Jij a b = ^{hi a hj a h* b h* b }M a Mb and Y^, is the Higgs abundance 
per isospin state. 

Sterile neutrinos also have lepton number, L = — 1 for left-handed N a fields 
and L = +1 for the conjugate N a . But what matters for particle propagation 
and transport of lepton number are the spin eigenstates rather than chiral states. 
When neutrinos are ultra relativistic, hclicity and chiral states are almost iden- 
tical, but when they are non-relativistic Majorana masses take over and lepton 
number vanishes. In fact, a positive (negative) helicity eigenstate has a well 
defined lepton number expectation value equal to the neutrino speed u (—it). 
Hclicity is not an invariant quantity, however, there is a privileged reference 
frame, the comoving thermal bath rest frame. In that frame, isotropy ensures 
that the spin density matrix is diagonal in the helicity basis. This is not true in 
any other frame, and means that neutrinos can be divided in two populations of 
opposite helicities. 

Under the non-equilibrium conditions of the leptogenesis era the two neutrino 
helicity populations develop unequal abundances and net lepton numbers in 
much the same way as the standard leptons. From the same decay, inverse decay 
and scattering processes one derives the Ni, lepton number source terms, 

(L b ) s = -r N ]T A^NJj - N b k) eq , (9) 

^-eq 

aij 

where tm is a coefficient estimated to be of the order but larger than unit. 
Summing both lepton and neutrino contributions, one obtains the total B — L 
asymmetry as a function of the temperature: 

where c « 2ci{tm — l)/7r ~ 0.1, Yf is the massless fermion thermal abundance 
and f(x) = 2 - (x 2 + 2x + 2) e - x . 

The other question is how fast Majorana masses erase lepton number. In a 
decay 4> — * l%N a , lepton number suffers a variation AL = 2 if N a has positive 
helicity, which occurs with a branching fraction of about 10 M 2 /T 2 . It makes 
an average variation AL w 20 M 2 /T 2 per decay and the N a production rates 
through Higgs decays are equal to HT' a /T with T' a w ^K a M a . Thus, lepton 
number violation enters in equilibrium at a temperature T in w ^ 20T' a M 2 /T 2 n , 
i.e., 

^3 10 
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The AL = 2 scattering processes mediated by Majorana neutrinos like U4> — ► lj<f>, 
that are so important when the Majorana neutrinos are non relativistic and 
even at lower temperatures when they are off-shell, are completely negligible 
in comparison with the Higgs boson decays AL — 2 transitions at relativistic 
temperatures. The relative strength is of the order of /i 2 a /200. For constants K a 
as large as 200, the sterile neutrinos enter in thermal equilibrium at temperatures 
T a w 0.35 K a M a ~ 70M a but the Majorana masses, i.e., the AL ^ transitions 
only enter in equilibrium at T- ln ~ 9 M$ ~ T a /8, where M3 denotes the largest of 
the Majorana masses. T; n is quite smaller than the relaxation temperatures T a 
and the larger the K a constants are the smaller is T; n in comparison with T a . 

4 Protecting B — L 

At temperatures above Tj n , B — L is not significantly damped and is slowly 
transferred via Yukawa interactions to the right-handed leptons and quarks, ej, 
Ui and di (with the exception of the top, ta, very early in chemical equilibrium). 
The dominant processes are Higgs boson decays into leptons ^1], 4> <-> eJi, 
kincmatically forbidden in the case of quarks due to comparable quark thermal 
masses ^5], and scatterings like 4>tR «-> bn<f), b^th <-» butu. These processes 
generate the otherwise zero asymmetries Q ei = L ei and Q Ui = B u . — B UR , 
Qdi = B^ — B UR . For right-handed quarks, the contrast asymmetries with re- 
spect to the up quark ur are the appropriate charges because, contrary to the 
individual flavor asymmetries, they are conserved by QCD instantons and 
are only violated by Yukawa couplings. Their growth rates are determined by 
the total Yukawa production rates, r x Yf, of x = e i, Ui — cr, di, and degeneracy 
parameters rj — \ijT as, 



neglecting flavor violation and ur Yukawa coupling. r x scale as the temperature 
and the ratios r x /H = T x /T define equilibrium temperatures, T x , under which 
the respective detailed balance equations are enforced. We estimate 2], ss 2 • 
10 12 GeV and T T w 10 12 GeV. On the other hand, weak sphalerons transfer part 
of the B — L asymmetry to the baryon sector and the baryon number growth 
rate |20l2ll^lll| in a volume V is given by 



where -T S ph is the sphaleron diffusion rate. 

The integration of the above equations requires knowing the relations be- 
tween the chemical potentials and the leptogenesis sources given in equations 
{7J and 10. In general there is no direct relation because the partial particle 



Qui 



2(77i s -V4> ~ Ve t )r ei Y f , 



(12a) 
(12b) 
(12c) 




(13) 
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asymmetries are not conserved due to flavor changing Yukawa couplings. As 
soon as the first neutrino, say N c , is in equilibrium (T < T c ), the heavy neu- 
trino and standard lepton flavors are rapidly violated enforcing their chemical 
potentials r\i i to be identical: % + ??jv c + ?70 = 0. The Higgs communicates also 
its asymmetry to tu and top quark doublet: r\ qt + 77^ — rjt R — 0. Because the Q 
asymmetries increase rapidly with time, initially as T , we concentrate on the 
latest period before Ti n , where sterile neutrinos are already in equilibrium but 
Majorana masses do not violate lepton number fast enough to yield null neutrino 
chemical potentials. Equations (|12a|l - 113|l can be written in the form 



where the coefficients a x are constant in a radiation dominated Universe but 
depend on the all set of chemical potential constraints. Assuming that Tin lies 
above 10 13 GeV one obtains for Qb and the baryon number, at ~ 0.48 and 
as ~ 0.24, respectively, while a T w because r\i i — r/^ — r\ ei ~ 0. The integration 
in time gives asymmetries per unity of entropy equal to 



where g(x) = 2 + (x + 4 + 6/x)e~ x — 6/x. 

The Q asymmetries grow with T -1 (g — > 2) down to the temperature T ln < 
T a where lepton number starts to be rapidly violated. Then, a new constraint is 
enforced, rji i + 77^ = 0, that would cause B — L and all the other asymmetries to 
vanish if there were no decoupled Q asymmetries. B — L falls down to the level 
of the dominant right-handed fermion asymmetries, or total baryon number, 
and is constrained to be a linear combination of them. They remain marginally 
damped as long as the respective Yukawa processes or weak sphaleron processes 
are slow, as implied by equations (|12a|) - (|13l) . As the temperature drops down a 
new chemical potential constraint is enforced whenever any of these interactions 
enters in equilibrium and B — L falls down further becoming a linear combination 
of the rcmanescent set of decoupled asymmetries. The final B — L asymmetry is 
determined by the Q asymmetries that are still approximately conserved when 
the B — L violating, AL — 2, processes become negligible. From then on B — L 
is exactly conserved. Assuming that this happens above 2 ■ 10 12 GeV, the 6r 
Yukawa and weak sphalerons are still out of equilibrium and the final B — L 
asymmetry is dominated by the Qb and B contributions: 



temperature T- m ~ 10 M3 and remain practically constant after that. The final 
B — L asymmetry is 



a x {B-L)-\Q 12 GeV, 



(14) 



dT- 1 




(15) 





(17) 
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As the temperature goes down and all right-handed leptons and quarks grad- 
ually enter in equilibrium B — L is redistributed among particles to eventually 
yield the present baryon number asymmetry Bo w ^(B — L) [13121)123] . One 
obtains B — L and baryon asymmetries at exactly the observed order of mag- 
nitude, B — L > 10 -10 per unit of entropy taking precisely the prompt decay 
ratio values indicated by the atmospheric neutrino evidence, K a = r a /H ~ 10 2 , 
without assuming any particular hierarchy between Yukawa couplings. 

In the solution described above B — L is mainly originated from the baryon 
number and asymmetry = Bb R — B UR . It assumes that AL — 2 reactions 
go out of equilibrium before the b^ and weak sphaleron processes enter in equi- 
librium at T « 2 • 10 12 GeV. As equation indicates this requires that 
the light neutrino mass scale m 2 is not larger than 0.2 eV 2 , which is still one 
order of magnitude higher than the atmospheric neutrino mass gap. If m 2 goes 
beyond 0.2 eV 2 , the charm cr replaces bn in the role of protecting B — L but 
Q* w 0.03 Ql is too small unless it is compensated by a strong hierarchy in the 
Yukawa coupling structure. 

5 Conclusions 

To conclude, leptogenesis can operate in a regime where all sterile neutrinos, 
including the lightest one, decay promptly when they vanish from the Universe. 
The light neutrino mass scale m 2 = ^ m 2 determines the natural order of mag- 
nitude of the final baryon asymmetry. First, m 2 values between the atmospheric 
neutrino mass gap and 0.2 eV 2 imply that the heavy neutrinos decay promptly 
with ratios K a — r a /H ~ 10 2 , and reach thermal equilibrium at temperatures 
10 times or more above their masses. During thermal production, B — L is gen- 
erated in both lepton and sterile neutrino sectors reaching the 10~ 7 level. Later, 
when the Majorana masses rapidly violate lepton number, B — L is effectively 
protected by the chemically decoupled bn quark and baryon number sector if, 
again, m 2 < 0.2 eV 2 , provided that the Majorana masses lie above the bn and 
weak sphaleron equilibrium temperature around 2 ■ 10 12 GeV. Then, B — L falls 
down 3 orders of magnitude, which meets exactly the observed baryon number 
asymmetry B without any further conditions. 

The dependence of B — L on the parameters of the theory is totally different 
from the delayed decay scenario. In the usual Mi <C M%, M3 hierarchical case, the 
lepton asymmetry generated in N\ decays vanishes in the limit of degenerate light 
neutrinos and |8I9| is essentially proportional to the atmospheric neutrino mass 
gap and the imaginary part ImfOfJ, where O31 is an element of the orthogonal 
complex matrix O = iv m^^U h M^ 1 / 2 that parametrizes the couplings h ia , 
and m v — U m U T is the diagonalized neutrino mass matrix. On the contrary, in 
the thermal production regime B — L does not vanish in the degenerate limit and 
goes with lm[(0^0) 2 3 ]. Equations ljl5|l and l|16(l yield bn and B — L asymmetries 
of the required 10~ 10 order for a wide range of parameters. They only depend 
on light and heavy neutrino mass ratios and matrix O. 
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